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(54) Clock offset estimation 

(57) A method is provided for estimating the clock 
offset between two clocks located at first and second sta- 
tions (A,B) respectively of a communications network. 
The method involves determining for each of a succes- 
sion of update events that provide data concerning the 
offset between the clocks: a first measure 6„ which is a 
noisy measure of clock offset; a second measure d n 
indicative of the power of the noise in 6 n ; and a third 
measure t n indicative of the time interval between the 
update event and a preceding such event. The measures 
derived from each successive update event are used to 



effect a Kalmanf iltering update process spanning Ksuch 
events, to provide an updated offset estimate. The 
Kalman filtering process may also be used to derive an 
estimate of the error V n of the clock offset estimate, the 
first time derivative of offset, and an error estimate for 
the latter. A level of noise is preferably injected into the 
Kalman filter to provide the filter with a tracking ability, 
the magnitude of this noise being determined on the 
assumption that convergence has been achieved for the 
offset estimate error 
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Description 

Field of the Invention 

s The present invention relates to a method of estimating the clock offset between two clocks located at first and 

second stations respectively of a communications network. 

Background of the Invention 

io The synchronisation of physically spaced clocks is required in a number of computer and communication applica- 
tions. Thus, for example, if a super-computer is to be built by interconnecting several computers in a network and then 
used to run applications that need a tight clock synchronisation, it is necessary to provide some mechanism for syn- 
chronising the clocks at the individual computers. Another example where clock synchronisation is required is where it 
is desired to analyse a sequence of events taking place across a network where each event has been time-stamped by 

is a different clock; such a situation might arise in relation to analysing a fault condition on a communication signalling 
network (for example, one operating according to the CCITT Signalling System No. 7 Recommendations). 

Protocols that perform clock distribution across computer networks have been available for a number of years. The 
current IAB standard is the Network Time Protocol, NTP (see, for example, D.L Mills : "Network Time Protocol (version 
3), Specification, Implementation and Analysis". Network Working Group, RFC 1305, March 1992). The NTP isageneral 

20 purpose clock distribution protocol that provides clock synchronisation to arbitrary computer networks and is actively 
used in the Internet. As such, its memory and processing requirements are low. For present purposes, the NTP can be 
considered as made up of a transmit process, a receive process, an update procedure and a local clock process. The 
update process generates the best offset available with reference to a selected reference clock by performing clock 
selection and combination. The local clock part uses the data generated by the update process to adjust the phase and 

25 frequency of the local oscillator with a PLL (Phase Locked Loop) which removes input data noise and provides a stable 
timescale between network updates. 

The NTP outputs three quantities: a clock offset with respect to a selected reference clock, an offset error (dispersion) 
with respect to the selected reference clock, and the current round-trip delay to this reference clock. The offset error or 
dispersion can be decomposed in turn into several components: 

30 

1 . - measurement errors due to inherent clock reading errors; 

2. - errors accumulated since the last offset update due to skew errors (here, "skew" refers to the first derivative of 
the offset with respect to time); 

3. - network errors due to delay jitter. 

35 

Normally, delay jitter errors dominate the dispersion budget. The NTP uses a type II PLL to minimise the delay jitter 
errors. However, PLLs have the problem that even if they work "weir in the sense of good tracking capabilities and 
acceptable errors, it is never known how much "better" one could do. 

It is an object of the present invention to provide an improved clock offset estimation method. 

40 

Summary of the Invention 

According to the present invention, there is provided a method for estimating the clock offset between two clocks 
located at first and second stations respectively, the method involving the steps of: 

45 

a) determining for each of a succession of update events that provide data concerning the offset between said 
clocks,: 

a first measure e n , for each update event, which is a noisy measure of said offset; and 
so - a second measure t n indicative of the relative time interval between the update event and a preceding such event; 

b) using the measures derived from each successive update event to effect a Kalman filtering update process, 
spanning Ksuch events, to provide an updated offset estimate, the Kalman filtering process being based on the 
following system equations: 

55 

X n + 1 -X n+ t n X' n ,n^0 
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where X n , X' n are respectively the clock offset and its first time derivative for the nth update event, and (H^) na0 
represents a noise function in the form of a sequence of independent random variables with known upper bound; 
the first time derivative of offset being represented in the Kalman filtering update process as: 

X 'n = Pn X n 

where: 

Pn = (Pa7 Pn,*) 

X n=( X n X n-K+l) 

p n k being determined for each latest update event (for example, by a least squares method) 

The system equations are not themselves in a form suitable for application of the Kalman recursions and it is the 
expression of X' n in terms of p„X n that permits application of the recursions. 

Preferably, the upper bound on the noise is determined for each updates event by deriving a third measure d n in 
respect of that event. 

The Kalman filtering process may also be used to derive an estimate of the error V n of the clock offset estimate, the 
first time derivative of offset, and an error estimate for the latter. Advantageously, a level of noise is injected into the 
Kalman filter such as to provide the filter with a tracking ability, the magnitude of this noise being determined on the 
assumption that convergence has been achieved for the offset estimate error V n 

In order to remove excessively noisy events from the process, the aforesaid third measure derived from each event 
is preferably compared with a threshold value and if this value is exceeded, the event is discarded. Advantageously, this 
threshold value is the average value of the third measure for at least the last M preceding update events where /Wis a 
positive integer. 

The update events used to collect data may take a variety of forms. In one arrangement, each update event, is a 
messaging event involving sending an outward message from the first station to the second station with the receipt of 
the outward message by the second station resulting in a corresponding return message being sent from the second 
station to the first. In this case, the method involves effecting a succession of messaging events, recording for these 
events the timings of message transmission and receipt at the two stations as judged against the corresponding said 
local clock, and determining from these timings for each messaging event: 

said first measure in dependence on the difference between the journey times for the outward and return messages 
associated with that event; and 

said second measure in dependence on the time interval between the messaging event and the preceding such 
event. 

The messaging events may be effected individually or in groups; in the latter case, for each group, all but the mes- 
saging event with the lowest value of the sum of outward and inward message journey times, are discarded. 

Preferably, the time between messaging events (or between groups of such events if the events are sent in groups), 
is set to a minimum value during the accumulation of timings on an initial K messaging events, and is thereafter deter- 
mined such that the error in offset estimate has grown with time only to a predetermined value before the next messaging 
event or group of such events, as applicable, is effected, provided that a certain minimum time has elapsed. 

Advantageously, for each new offset estimate, a check is made as to whether the change in estimate between the 
new estimate and the preceding one is within a predetermined limit, the inter-event or inter-group time, as applicable, 
being set to a minimum if this limit is exceeded. 

Brief Description of the Drawings 

A clock-offset estimation method according to the invention will now be described, by way of non-limiting example, 
with reference to the accompanying diagrammatic drawings, in which: 

. Figure 1 is a diagram showing two stations A and B inter-communicating via a network, each station having its 
own clock; 

. Figure 2 is a diagram illustrating the timing of a message event as seen by the clocks at the two stations A and 
B;and 

. Figure 3 is a flow chart illustrating the control effected over the time interval between the sending of bursts of 
messaging events. 
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Best Mode of Carrying Out the Invention 

Figure 1 illustrates the generalised situation to which the present clock-offset estimation method is applicable. Two 
stations 1 0, herein referred to as stations A and B, inter-communicate over a network 1 1 . By way of example, the stations 
5 A and B may be computer workstations and the network 1 1 a LAN (Local Area Network) or a WAN (Wide Area Network). 
Each station A, B has a local clock 12, herein designated clock A for station A and clock B for station B. In order to 
enable a correlation of events at the two stations A, B, it is necessary to know the offset of one clock 12 relative to the 
other. This offset will generally change with time and the first time derivative of the offset will hereinafter be referred to 
as the "skew"; the second time derivative of the offset will be referred to as "drift". 
w The clock-offset estimation method to be described enables an estimation to be made of the offset and skew between 
clocks A and B. Generally, one station (in this case, station A) actively carries out an estimate of the offset of its clock 
A relative to the clock 12 (in this case, clock B) of the other station (station B); station B may be considered as being a 
"time server" and station A as being a client of B. 

This estimation process involves station A effecting a succession of messaging events between itself and station 
is B. Each messaging event involves the sending of an outward message from station A (source station) to station B 
(destination station) and the return of a return message from station B to station A. Timings associated with at least 
some of these messaging events are recorded and then processed to derive a clock offset estimate. 

In functional terms, station A can be considered as including, in addition to clock A, five blocks for effecting the 
estimation process. The blocks are a message event block 13 for sending and receiving message-event messages, a 
20 network load monitor 1 4 for discarding "noisy" messaging events, an offset and skew estimation block 1 5, a jump detec- 
tion block 1 6, and a coordination block 1 7 for coordinating the operation of the other blocks. Where station A is a computer 
workstation, these blocks wiirgenerally be implemented by appropriate processes run on the workstation processor, the 
exception to this being the network interface hardware that conceptually forms part of block 13. 
The functional blocks 13 to 17 will now be described in turn. 

25 

Messaging Event Block 13 

The messaging event block 1 3 is arranged to effect a series of messaging events, each involving sending a message 
and receiving the return message generated by station B. The time of sending of the outward message and of receipt 

30 of the return message are noted according to clock A whereas the time of receipt of the outward message and of sending 
back of the return message are noted according to clock B. General mechanisms for providing such a facility are known 
in existing network protocols such as the echoing and time-stamping mechanism of the ICMP protocol associated with 
the TCP/IP protocol suite; with such a mechanism, an outward message is time-stamped on sending by-the clock at the 
sending station (station A) and the message is then echoed back, as a return message, from a remote station (station 

35 B) which appends its own time-stamps corresponding to the time of receipt and of echo back at the remote station. On 
receipt of the echoed-back message, the original sending station has the required timing information in the received 
message. It will, of course, be appreciated that the timing information can be passed from each station A, B to the means 
effecting offset estimation processing (whether or not part of station A) by mechanisms other than including time stamps 
in the outward and return messages. 

40 In the present example, it is assumed that the stations A, B run the ICMP protocol and use its mechanisms to effect 
and time the messaging events, station B providing timestamps for both message receipt and echoback. 

Referring now to Figure 2, the flow of a messaging event (taken as the n th such event) is illustrated relative to the 
time lines 20, 21 for clocks A and B respectively. Thus, an outward message (represented by line 22) is first sent from 
station A to station B which then sends back a return message (line 23) to station A. The time of sending of the outward 

45 message 22 is recorded as being at time T n by clock A (the suffix "n" indicating that the timing relates to an n th messaging 
event); note that in practice, time T n is slightly before the outward message is actually transmitted from station A, this 
interval being a handling interval needed by station A to effect transmission after timestamping at T n . The receipt of the 
outward message 22 at station B is recorded by a time stamp at time R n according to clock B (note that for clarity, timings 
according to clock B are primed); between receipt of the outward message 1 7 and R n is a small handling interval. The 

so sending by station B of the return message 23 is recorded by a time stamp at time S' n ; between S' n and transmission 
of the return message 23 is another small handling interval. The time of receipt of the return message at station A is 
recorded as time A n according to clock A, there again being a small handling interval between receipt of the return 
message and time A n . The messaging event block 13 thus has available time values T n , R' n , S n , A n following the n l 
messaging event (note that if station B only provides one time stamp for message receipt and echoback, then R' n = S' p ) . 

55 Generally, of course, there will be an offset between clocks A and B. Thus, when clock B records times R' n , S' n , 
clock A will show different values R n , S n respectively. More particularly, if X(t) is the clock offset at source time t (clock 
A time), then: 

R n =R' n +X(R n ) (D 
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S n =S' n +X(S n ) (2) 

Similarly, times A n and T n at the source (station A) correspond to times A' n and T' n at the destination (station B). 
For the n th messaging event, block 13 calculates the quantities: 

5 ' . 

(A n - S'J+ (R' n - T n ) 
Half round trip delay d n = — n/ 1 n ^ (3) 

"> Packet Offset (A ° ' SV " ffl ' n ' Tn) (4) 



The quantities d n and 6 n are be directly computed from the timestamps. 

The larger d nt the "noisier" the offset information of the n th messaging event In order to reduce the noise (that is, 
is to reduce the half round trip delay d n ), rather than sending isolated messaging events, block 13 is arranged to send 

bursts each of several messaging events, with, for instance, the events in a burst being separated by about one second. 

The first messaging-event packet in a burst acts as a "warning shot" and has the effect of moving the destination address 

into the cache of routers in the network. As a result, subsequent messaging-event packets of the same burst normally 

end up with a smaller half round trip delay. For every burst of messaging events, block 13 selects the messaging event 
20 with the smallest value of d^ 

Thus, for the n th burst of messaging events, the block 1 3 selects the messaging event with the smallest value of d n 

and treats this selected messaging event as the n th messaging event; block 13 then proceeds to output the following 

values for this selected event: 

25 d n Q n T n A n 
Network Load Monitor 14 

The network load monitor 14 tracks the network load. The network load is estimated from the half round trip delay 
30 values. Upon output of the n th messaging event data by block 1 3, the load monitor 1 4 first uses the half round trip delay 
value d n to update an admission threshold A to a new value A n . Preferably, A n is computed as the average half round 
trip delay seen over the last M selected messaging events from the last M messaging-event bursts. It is the value of M 
rather than which average (e.g. arithmetic mean, median,...) that plays an important in monitoring the network load, a 
large value of M making the monitor less adaptable and a small one making the monitor redundant (M - 1 would mean 
35 that each messaging event would always be accepted). 

After the threshold value A n has been updated, the value of d n is compared to the threshold. If d n is less than A n 
then the n th selected messaging event is considered valid and the values output from block 13 are passed to the esti- 
mation block 1 5 for clock offset and skew estimation ; otherwise, the n th messaging event is discarded and the coordination 
block 17 informed. In this latter case, the coordination block 1 7 initiates a new messaging event burst by the block 13. 

40 

Offset and Skew Estimation Block 15 

Block 1 5 estimates the clock offset X(t) and the clock skew X'(t) from the valid selected messaging event observations 
(d 0 , 6 n , T n , A^, nsO. This is carried out using Kalman filtering as described below. 

45 

Clock Offset Model. Assumptions and Notation 
Evolution Equation: 

so It is assumed that the clock offset X(t) at source time t has an evolution equation of the form: 

~\ 

t 

X(t) = x+ ^X'(s)ds (5) 
b 

55 

the quantity X(0) being called the initial clock offset and the function X'() being the clock skew (first time derivative of 
offset). 
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Assumptions: 

(1) - the drift (that is, the second time derivative of offset) is practically zero so that between two consecutive valid 
messaging events, the clock skew remains practically constant; 
5 (2) - the clock skew per unit time is assumed to be very small (a few parts per million): 

\X'(t)\ «1. f a 0 

(3) - protocol stack processing and routing delays (source-destination; destination-source) are on average symmet- 
70 rical. 

Assumption (3) is often not met in practice; however, this is not critical as this only leads to a constantly biased offset 
value, the bias being dependent on the degree of asymmetry. In contrast, assumptions (1) and (2) are more critical in 
nature (but see later for how a failure of assumption (1) can be handled). 

15 ' 

Notation: 

Let: 

20 X n = X((A n +T n )/2) 

X' n = X'(Q, (A n +T n )l2 sf< (A n + 1 +T n + 1 )t2 

Basic Equations 



25 



Single-Event Equation - The first equation that can be set down for relating messaging event data to the clock 
offset model, concerns each event taken in isolation. The quantity 6 n (which is half of the difference between outward 
and inward packet journey times) would actually be a true measure of clock offset were it not for the variations from the 
symmetry assumptions noted above that affect any single messaging event. It may therefore be said that: 



30 



e n oX„ + W n ,nsO (6) 



where {W n ) n ± 0 represents a sequence of independent variables of zero mean (if assumption 3 holds) bounded by d n 
and with variance (a„) 2 . W n is called the measurement noise. The initial value X 0 is a random variable of mean 8 0 and 
35 variance (co) 2 . 

Event-Sequence Equation - The second equation that can be derived for linking the measurement results with the 
offset model, concerns changes taking place from one messaging event to the next. Clearly, the offset at event (n+1) 
will be related to the offset at event (n) by the clock skew (first time derivative of offset); however, there will also be some 
randomising factor involved. The second equation has the form: 

40 

X n + i = X n+ t n X- n+ W n , v ^0 (7) 

where: 

45 t n =(A n + 1 +T n + 1 )/2-(A n +T n )I2,n*0 (8) 

and ( W n )n2i represents a sequence of independent zero mean random variables with variance (v n ) 2 . W' n is called the 
driving noise. The random variable sequences W n and W' n are independent of each other. 

so Kalman Filtering 

Based on equations (6) and (7) above, a Kalman filter is provided to recursively compute the clock offset estimate: 

X n = E[X n |Oo,...,0n], n>0, 

55 
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and the mean square error: 
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V n 2 = E[(X n - I) 2 ], n>0. 



In this process, the clock skew (the first time derivative of offset) is estimated from the measurement data. More generally, 
there may be situations where it is important to capture the offset dependency on the skew, the drift and higher derivatives 
10 (this would permit assumption (1) above to be breached); in the following it is shown how this can be done with a K 
dimensional Kalman filter though, for clarity, it is assumed that only the skew is needed (the same approach being usable 
for higher order time derivatives). 



15 



Offset Estimate - First, the skew is written as: 



X 'n = Pn X n 



(9) 



20 



where: 



Pn = (Pn.r P/) f tf) 

X n = ( X n X n-K+l) 



(10) 
(11) 



The system equations (6) and (7) can now be rewritten in the following form (superscript T denoting the transposition 
25 operator): 



30 



where: 



e n =H"X n+ W n . neO 



(12) 
(13) 



H= (1,0 0) 



(14) 



35 



40 



( 1 + a n> i a n , 2 
1 0 



F n = 



0 

0 
0 



1 

0 
0 



<*n,A'-l °tn,K \ 
0 0 



0 

0 



0 
0 

o / 



(15) 
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50 



If 



W' n+n =(W' n + 1 ,0 0) T 

a a/ = 'nPn./ 



(16) 
(17) 



denotes the estimate of X n given 6 o 6^, and the matrix C n denotes the Kx Kcovariance matrix of X n , that is: 

55 then the theory of Kalman filtering provides formulas for the updates of C n and X n at time v 
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These updates are as follows: 



+ / (w, T + a,*,)*/ \ 

1(<W> 2 +("„ + ;) 2 + V„ 2 + a|1 2 J 



(19) 



' tP n F .T n . {F n C n F n J + Q n+1 )HH T {F n C n F n + Q n+1 ) 



where: 



V n 2 = H J C n H (21) 

a n 2 = a n C n a n T (22) 

0 n-( 0 ft»-- a n,#f) ( 23 > 

Q n+ ,=diag{^ 2 ,0 0} (24) 

The initial values are set after the first K messaging events as: 

A 

•^o = (Q. k+l ,....,Q 0 ) T (25) 



(20) 



C 0 = diag{(a. K+7 ) 2 (a/} ' (26) 

Q 0 = diag{(y 0 ) 2 ,0 0} (27) 

The offset- estimate at time (A n + , + T n + , J/2 is: 

X n+1 = X n+ t n X- n + {J n+1 ) 2 (e n+1 -{X n+ t n X' n ))nzO (28) 
and its mean square error is: 

(V n+1 ) 2 = (c n+1 ) 2 (J n+1 )\n^0 _ ' (29) 

where: 

Skew Estimation - To enable the clock offset to be derived according to the above formulae, the vector p n is estimated 
using a least squares estimate. If E = (e,y) is the {K-1)xK matrix given by : 

1if/=y (31) 
ej = -1 if ;'= ■•+ " 
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0 otherwise 

where 1 ^&(K- 1) and 1 then the (K-1) vector with i-th component, 

Is/sKcan be written as T n X n where T n = diag{1/f n _, Mt n . K ^ 1 }E n . 

The covariance matrix G n of the random vector T n X n is equal to: 

G n =T„C n 7-„ r (32) 

The function to minimise in order to obtain a least squares estimate of p and of X n is: 

X 2 = (T n X n -X' n e?G n '(T H X n -X' n e) (33) 



where G„ # is the generalised inverse matrix of G n and e = (1, ....1) 7 ! From this it can be shown that the clock skew 
20 estimate is then: 



25 




(34) 



30 



and thus, 



e T G„V„ 



a n ' n 

e G n e 



(35) 



35 



and the skew variance is given by: 



40 



(a*.) 2 =Var(X' n ) = ^ # 



(36) 



Dynamic Behaviour - There is a trade-off between the Kalman filter tracking ability and the Kalman filter error 
estimate. The driving noise power {v n ) 2 is used to control this trade-off. If the driving noise is too large, V n 2 tends to a 
value that is unneeessarily large. The question to address is what is the right power of driving noise that needs to be 
45 injected into the Kalman filter to preserve its tracking ability without creating error estimates that are unnecessarily large. 

The injection noise level used is that which would correspond to having achieved error convergence, that is, 
v n = v n+i • Then, at the reception of the (n+1)th messaging event, the level of noise to be injected into the Kalman 
filter is given by: 



50 



(37) 



55 where = max {0,x} . 

This approach seems to enhance the tracking capability of the filter at the expense of the offset error performance. 
This is justified since the errors obtained by the Kalman filter are normally an order of magnitude below the half-round 
trip delay 
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Implementation 

For the (n+1)th valid messaging event, the estimation block 15 has the following new data to work with: 
d^t - the half round trip time which provides an upper bound on the measurement noise; it is generally used to 
5 determine the standard deviation of the measurement noise when evaluating the above formulae; 
e^t - the packet offset; 

t n - the inter-event time (derived according to equation 8 above). 
This data is then used to update the offset estimate as follows: 

w 1) the matrix F n is determined by: 

- calculating the vector p n using formula (35); 
calculating a a/ using formula (1 7); 

- calculating Pausing formula (15); 



15 



40 



45 



2) the matrix is determined by: 



- calculating the driving noise power ( d 2 using formula (37) , the value of V n having been determined in relation 
to the previous messaging event; 

20 - calculating Q m / using formula (24); 

3) the value of a* is determined using formula (22), the value of C n having been determined in relation to the 
previous messaging event; 

4) the offset is estimated using formula (19); 

25 5) the mean square error of the offset is estimated by: 

- determining the value of (J^?) 2 using formula (30); 
determining the mean square error using formula (29); 

30 6) a skew estimate is derived using formula (9) together with an estimate of its mean square error using formula (36); 
7) the value of 1 is determined using formula (20). 

It will be appreciated that the numerical computations of the K dimensional Kalman filter need to be done carefully 
in order to avoid numerical problems (such as division by zero); appropriate techniques for avoiding such problems will 
35 be apparent to persons skilled in the art. 

Jump Test Block 16 

A highly unexpected value of e n may indicate a drastic change in the network behaviour or conditions; in particular, 
it may indicate a jump in the value of clock offset. Clock offset jumps will generally be rare and much larger than fluctu- 
ations due to noisy data. Block 1 6 tests for such jumps and if one is detected, block 1 6jnforms the coordination block 
1 6 which sets the inter-burst time to minimum to obtain further data as quickly as possible. 

The jump test carried out by block 16 also has the function of ensuring that the Kalman filter maintains its tracking 
abilities - if it does not, then this will be picked up by the jump test. 

The nature of the jump test is as follows (this being one example only, other tests being possible). On the basis that 
the clock skew will only vary slowly with time, the hypothesis: 

H 0 :X n+7 = X n W n X' n (38) 

so against the hypothesis: 

If hypothesis H 0 holds, then it is taken that no offset jump has occurred. Under hypothesis H 0 , it may be expected that: 

e n+r «HV n *„+W'n + t)+Wn + r (40) 



55 
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and since e^. ? 

r n +, = 8„+ ; - H T F a x a (4i) 
>/(*„ + ,) 2 +("„ + ,) 2 + K 2 + «„ 2 



and X n are independent, it follows that the quantity: has zero mean and unit variance; hypothesis Ho can be rejected if, 
for example, \ Y n+1 \ > z, where z is a pre-specified value. 

The numerator in formula (41) is the absolute size of the clock offset jump. Thus if one is interested in detecting 
jumps of a fixed size, regardless of the normalisation introduced by the denominator of formula (41 ), it suffices to compare 
the numerator in formula (41 ) with the fixed size threshold of interest. 

If hypothesis H 0 is accepted no action is taken. However, if hypothesis H 0 is rejected, the coordination block 17 is 
informed. 

Coordination Block 17 

The coordination block 1 7 coordinates the operation of the other blocks and, in particular, controls the time between 
message-event bursts in dependence on how well the estimation process is doing according to feedback provided by 
the other blocks. Clearly, at start up when data on K valid selected messaging events is required to start off the Kalman 
filter, the inter-burst time is preferably set at a minimum; however, once a steady state condition has been achieved, 
then the inter-burst time will depend on the allowable level of error for the offset estimate (the longer the time between 
bursts, the greater the potential error). More particularly, the offset error t time units after time {A n + T n )l2 is given by: 

V 2 (t) = (H T +tp n )C n (H T +tp n ) T (42) 

where: 

V 2 (0)= V n 2 = H T C n H (43) 

is the offset error at time ( A n + T n )/2 

The offset error increase will typically be controlled not to exceed a value ed n where e is a preassigned quantity. If 
(e d n - V n ) > 0 , the (n + 1 )th burst will be sent t time units after the nth burst, where t is the positive solution of: 

(a k .) 2 t 2 + 2bt-{e 2 d n *-V n 2 ) (44) 

^ n 

where 

b = fi n C n H=H T C n p n T 

Since usually 

(o k .) <<b 

and since also there will be a minimum allowable time t min between bursts, the inter-burst time f will be set such that: 
- if (z 2 d n 2 - V n 2 ) >0, f is the greater of: 

/ 2 2 2. 
pC n H + H T C n p n T 

and 

otherwise, t = t min 



Figure 3 is a flow chart illustrating the operation of the coordination block 1 7 in controlling the sending of message-event 
bursts by block 13. At startup, a counter has its value c set to zero (step 30); this counter is used to count the number 
of accumulated valid messaging events (Ksuch events being required to start off the Kalman filter). Next, the inter-bust 
time f is set to f min (step 31) and after time f (step 32) the coordination block 1 7 triggers the block 13 (step 33) to send 
5 a messaging event burst. 

If block 13 receives no return message from station B then this is reported to the coordination block 17 and block 
' 17 returns to step 31 (see step 34). 

If however, at least one return message is received by block 1 3 for the messaging-event burst, block 17 is informed 
that a messaging event has been effected and block 1 7 then proceeds to check the count value c (step 35) If c is less 
w than the required number Kneeded to run the Kalman filter, the value of c is incremented (step 36) and block 1 7 returns 
to step 32; however, if c equals K, then block 1 7 awaits a "valid event" indication from the network load monitor 14. 

If the network load monitor 14 reports that no valid messaging event results from a burst (tested in step 37) . then 
block 17 returns to step 31. However, if a valid messaging event is reported, block 17 next awaits receipt of the new 
' offset, skew and error estimates from the estimation block 1 5 (step 38). 
is Thereafter, if block 16 reports that no offset jump has been detected (tested in step 39). block 1 5 proceeds to output 
the offset estimate (step 40) and determine a new value for f (step 41) before returning to step 32. In the presence of 
an offset jump, block 15 returns to step 31. 

It should be noted that the offset estimate determined by block 15 is the offset at time (A n +T n )/2 . The offset 
estimate at time A, is determined by taking account of the offset change due to the value of the skew over the period 
20 {A n - T n )/2 . The offset at any intermediate time between two estimates can be similarly derived. 

Variance Estimation 

In the foregoing, the variance ct„2 of the measurement noise W„ was equated to its known upper bound d„2. However, 
25 thequantity introduced in formula (41) must have a variance of 1 which is not achieved with a n+ ,»d n+ , because 
o , is too high. It is. in fact, possible to define a parameter y, between 1 and 0 in value, such that if a n+ , = yd n + , then 
the variance of V^, is equal to 1. Furthermore, in this case, all variances are scaled by y; as a result, the offset error 
V/ which was previously derived using a „ 2 =d„ 2 , now becomes scaled down to y 2 ^ 2 . 

The value of y is estimated for the nth valid messaging event as the running variance of all the quantities V„ to date. 
30 where the quantities ( vg comprise the sequence of values obtained from formula (41 ) with o n = d n . 



Variants 



It will be appreciated that many variants are possible on the foregoing. For example, the steps taken to exclude 

35 noisy messaging events may be omitted, though again this is not preferred. 

Although in the described embodiment, messaging events have been used to collect the measurement data to which 
the Kalman filtering process is applied, it is possible to collect the necessary update data in other ways. For example; 
one station might periodically inform the other of what it thinks its offset is; alternatively, in appropriate circumstances 
one could infer the offset of a station by the way it generated traffic. The update measures required are an offset measure 

40 e„ (which will inherently be noisy), an upper bound d n on the noise of the offset measure, and a knowledge of the time 
between successive update measures. With regard to the upper bound on the noise power in the offset measure, if the 
environment concerned is sufficiently well known, it is feasible to dispense with making an upper bound measure for 
each update and simply use a predetermined value; however, this is not preferred. 

45 Claims 

1 . A method of estimating the clock offset between two clocks located at first and second stations respectively of a 
communications network said method comprising the steps of: 

so a) determining for each of a succession of update events that provide data concerning the offset between said 

clocks,: 

- a first measure for each update event, which is a noisy measure of said offset; and 

- a second measure t n indicative of the relative time interval between the update event and a preceding such 
55 event; 

b) using the measures derived from each successive update event to effect a Kalman filtering update process, 
spanning K such events, to provide an updated offset estimate, the Kalman filtering process being based on 
the following system equations: 
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= W n .n*0 

s where X„, X'„ are respectively the clock offset and its first time derivative for the nth update event, and ( WJ nj:0 

represents a noise function in the form of a sequence of independent random variables with known upper bound ; 
the first time derivative of offset being represented in the Kalman filtering update process as: 

10 

where: 

Pn = (Pn,1 .PnX> 

15 Xn = 0<n *n4C+l) 

p n ^ being determined for each latest update event. 

2. A method according to claim 1 , wherein said upper bound on the noise is determined for each update event by 
20 deriving a third measure d n in respect of that event. 

3. A method according to claim 2, wherein said Kalman filtering process is also used to derive an estimate of the error 
V n of the clock offset estimate. 

25 4. A method according to claim 3, wherein a level of noise is injected into the Kalman filter such as to provide the filter 
with a tracking ability, the magnitude of this noise being determined on the assumption that convergence has been 
achieved for the offset estimate error V n 

5. A method according to claim 2, wherein an estimate is also made of the value of the first time derivative of the offset 
30 by using covariance data obtained from the Kalman filtering process. 

6. A method according to claim 1 , wherein between offset estimates, the actual current clock offset is estimated by 
combining the latest offset estimate with a component dependent on the lapsed time since that estimate and the 
first time derivative of the offset. 

35 

7. A method according to claim 2, wherein a said update event is discarded if the said third measure derived therefrom 
exceeds a predetermined value. 

8. A method according to claim 7, wherein said predetermined value is the average value of said third measure for at 
40 least the last M preceding update events where M is a positive integer. 

9. A method according to any one of the preceding claims, wherein each said update event is a messaging event 
involving sending an outward message from said first station to said second station with the receipt of said outward 
message by the second station resulting in a corresponding return message being sent from the second station to 

45 the first; the method involving effecting a succession of said messaging events, recording for these events the timings 
of message transmission and receipt at the two stations as judged against the corresponding said local clock, and 
determining from said timings for each messaging event: 

said first measure in dependence on the difference between the journey times for the outward and return mes- 
so sages associated with that event: and 

said second measure in dependence on the time interval between the messaging event and the preceding such 
event. 

10. A method according to claim 9, wherein said messaging events are effected in groups each of at least two messaging 
55 events, and for each group, all but the messaging event with the lowest value of the sum of outward and inward 

message journey times, are discarded. 

1 1 . A method according to claim 9 or claim 10, wherein the time between messaging events or between groups of such 
events if the events are sent in groups, is set to a minimum value during the accumulation of timings on an initial K 
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messaging events, and is thereafter determined 
such that the error in offset estimate has grown with 
time only to a predetermined value before, the next 
messaging event or group of such events, as appli- 
cable, is effected, provided that a certain minimum 
time has elapsed. 

12. A method according to claim 11, wherein for each 
new offset estimate; a check is made as to whether 
the change in estimate between the new estimate 
and the preceding one is within a predetermined 
limit, the inter-event or inter-group time, as applica- 
ble, being set to a minimum if said limit is exceeded. 
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